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Abstract 

Within the framework of nonrelativisitic QCD factorization formalism, we present the next-to- 
leading-order relativistic corrections to T exclusive decay into gc plus J/b- The double charmonia 
can be prodnced through several immediate channels, i.e., T — g*g*g* gc + J/f’-, T ^ g*g*g* —>■ 
gc + J/'f’i and T —)■ 7* —>■ r/c + J/ij- The amplitudes of these three channels are obtained accurate 
up to 0{aa1v'^), and 0{aagv‘^), respectively, where v indicates the typical heavy quark 

velocity in bottomonium and/or charmonium rest frame. The decay rates are also presented. 
We find that the next-to-leading-order relativistic corrections to the short-distance coefficients 
as well as the decay rates are both significant and negative, especially for the corrections from 
bottomonium. More seriously, the decay rates are even brought into negative by including the 
relativistic corrections, which indicates the poor convergence for the velocity expansion in this 
kind of process. Detailed analysis is given in the paper. 
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I. INTRODUCTION 


Nonrelativistic QCD (NRQCD) [1] is a powerful and successful effective field theory in 
describing the quarkonium production and decay. The processes of bottomonium exclusive 
decay into double charmonia are interesting channels to investigate quarkonium, which have 
been extensively studied both in theory and experiment [2, 3]. The r]h decay into double 
J/ijj was studied in Refs. [4-8], and the radiative and relativistic corrections to this process 
were investigated subsequently in [9-11]. The P-wave bottomonium exclusive decay into 
double charmonia was also studied in Refs. [12, 13]. The relativistic corrections to these 
processes were obtained in Refs. [11, 14], and the radiative corrections to xm J/'’PJI'’P 
were calculated recently [15]. For the case of T, the T —)■ rjcJ/pJ and T —)■ J/ipXcj were 
separately studied in Ref. [16] and Ref. [17]. 

Comparing with other bottomonium exclusive decay processes mentioned above, T —)■ 
rjcJ/'P has some distinctive features. First, there is no tree-level pure QCD Feynman dia¬ 
gram. Second, this process is further suppressed by a helicity selection rule due to nib ^c- 
As a consequence, one may expect that the decay rate is relatively small [16]. Nevertheless, 
it is still meaningful to investigate this process based on the following reasons. First, there is 
little study of the relativistic corrections to the quarkonium decay and production processes 
with the vanishing tree-level Feynman diagram, so it is theoretically interesting to study 
the relativistic corrections and the convergence of the velocity expansion for these processes. 
Second, it has been known for a long time that the relativistic corrections to J/p) LH 
and J/p) —)■ 37 are both signihcant and negative [18], which even surpasses the contribu¬ 
tions at leading order. The study in Ref. [19] indicated that the 0{asV^) corrections to 
J/p) ^ ?)X are also huge; meanwhile, the study in Refs. [20, 21] indicates that the relativistic 
corrections to T —)■ g*gg —)■ ccgg are extraordinarily large, despite the small 0{v'^) NRQCD 
matrix element for T. The common feature for these processes is that T hrst annihilates 
into three gluons (photons). For the process T —)■ rjcJ/P’t the main contributions come from 
the immediate channel T —)■ g*g*g* —)■ rjcJ/P’ [16]. As a consequence, it is compelling to 
know whether large relativistic corrections are also appearing in T —)■ g*g*g* . Third, there 
has been a longstanding “p-7r” puzzle in quarkonium decay [22, 23]. The structure of Feyn¬ 
man diagrams for p) ^ pn is rather similar to T —)■ rjcJ/P^i so one may wonder whether 
the puzzle originates from the relativistic corrections for J/p). Based on these reasons, we 
will investigate the relativistic corrections to T —)■ rjcJ/P’ from both the charmonia and the 
bottomonium. 

Up to the next-to-leading order in velocity expansion, there are three channels which 
contribute, T —)■ g*g*g* — )■ rjcJ/P’, T — t g*g*'y* —t PcJ/i’-, and T —)■ 7* —)■ rjcJ/P’- The 
leading-order (LO) amplitudes for the three channels in powers of o;* and a scale as o;^, aa^, 
and aas, respectively. Though the second channel is suppressed by a factor a/as relative 
to the hrst channel, there is a fragmentation enhancement (7* —)■ J/p>) in T —)■ g*g*Y 
TjcJ/P^i as a result, we also consider this channel with J/pj produced by a photon in our 
calculation. Moreover, since T ^ 7* —)■ PcJ/P’ is of a similar Feynman diagram structure 
as rjcJ/P^i th® radiative corrections to which were found to be large [24, 25], the 

calculation in this channel will be accurate up to O^aalv"^). ^ 


^ We do not consider the pure QED contributions in T —>• 7 * —>• gcJ/P, due to the fact the pure QED 
contributions are less than a third of other contributions [26], and what is more, T —> 7 * —>■ gcJ/p 
contributes merely a small fraction in T —>■ gcJ/P decay [16]. 
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The rest of the paper is organized as follows. In Sec. II, we describe the NRQCD formula 
for T exclusive decay into double charmonia and present the kinematic descriptions and 
required techniques. The short-distance coefficients (SDCs) are obtained in detail in Sec. III. 
Section IV is devoted to numerical predictions. Discussions and summary are also presented 
in this section. 


II. NRQCD FORMULA AND KINEMATIC DESCRIPTIONS 

According to Lorentz invariance, we are able to factorize T exclusive decay into double 
A-wave charmonia as 


A = (r/,J/^(e2)|T(ei)) = (1) 

where e tensor indicates the Levi-Civita symbol, ei and 62 represent the polarization vectors 
of the initial state T and final state J/'^, Pi, and P 2 denote the momenta of r]c and J/'^, 
respectively. 

As mentioned in Ref. [11], the exclusive decay of a bottomonium into double charmonia 
involves the annihilation of a bb pair followed by the creation of two pairs of cc. One may 
guess that the generalization of the NRQCD factorization for the electromagnetic decay 
or light-hadronic decay into this exclusive mode is possible. So according to the NRQCD 
factorization formalism, the Lorentz scalar A can be factorized as 

A = ( 0 ),.( 0 )j/*{C>)t Co (l + + C2,2(w^),. + C2,3{v'^)j/t + 0 {v‘)) , 

( 2 ) 

where Cq and C 2 correspond to the SDCs, which can be calculated perturbatively in power 
of as- For simplicity, we dehne shortcuts for the matrix elements 


(0}t = |{0|xV.e;V.|T(£.))p, 

(3a) 

{0)„ = 

(3b) 

{0}j/,l, = |{J/'i;i(£2)|V’'f ■ £2X|0>t. 

(3c) 

(0|V<T.£:(-ih’)V|T(£l)> 

(4a) 


(4b) 

. , 2 \ _ (W(e2)|^^cr-e2(-|^)\|0) 

ml{J CT ■ e2x\0) 

(4c) 


In Eq. (2), the bottomonium is normalized nonrelativistically on both sides; in contrast, the 
charmonia are normalized relativistically on the left-hand side, but nonrelativistically on the 
right-hand side. The explicit factor is responsible for the discrepancy. 

The SDCs are free of nonperturbative effects and can be obtained through the matching 
technique. To calculate the SDCs, we are allowed to substitute the quarkonium states to free 
heavy quark pairs of the same quantum numbers as the corresponding hadrons. Now A and 
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the matrix elements in ( 2 ) can be calculated perturbatively; therefore, the SDCs are readily 
derived. The spin-singlet and -triplet states of a free heavy quark pair can be extracted by 
employing the spin projectors in Ref. [27]: 


n3(Po,Po) 

ni(pi,pi) 


n 3 (P 2 ,P 2 ) 


Sy/W\,E{qQY[E{qQ) -h rm,] 

-i-( 

A^,E{q^)[E{qr)+m,y 

1 

Ay/2NlE{q2)[E{q2) + m. 


{j)o + mb)[pQ + 2E{qQ)]ii{fQ - mb), 
pi - me) 'y5[pi + 2E{qi)]{pi + mc), 
-{p 2 - me) f 2 [p 2 + 2 E{q 2 )]{p 2 + me), 


(5a) 

(5b) 

(5c) 


where a color-singlet factor is included, and E{qi) = with m = m?, for i = 0 

and m = me for i = 1,2. In (5), we take P* and to be the total momenta and half of the 
relative momenta of the quark pairs in the quarkonia, where i = 0 denotes T, i = 1 denotes 
rje, and i = 2 denotes J/y. Therefore, the momenta of the corresponding quarks pi and 
antiquarks pi in the hadrons are expressed as 


Pi = ^Pi + Qi, 
Pi 2^* 


( 6 a) 

( 6 b) 


Pi and qi are chosen to be orthogonal: Pi ■ qi = 0 . 

To obtain the next-to-leading-order (NLO) relativistic corrections for S'-wave quarkonia, 
we can expand the amplitude in powers of qi, keep the quadratic terms and make the 
following substitution 


LL 1/ 

QiQi 


^ I - 17^" + 


pi^pu 


AE{qP)Jd-l 


Qt 


( 7 ) 


where d is the dimension of spacetime. In addition, the contributions from E{qi) can be 
readily derived by expanding 


E{qi) = Jm? + q 


= m + + 0{^^) 

2m m^ 


( 8 ) 


where m = mb for i = 0 and m = me for i = 1 , 2 . 


III. THE RELATIVISTIC CORRECTIONS TO THE SDCS 

In this section, we calculate the SDCs of T —)■ peJ/y accurate up to NLO relativistic 
corrections following the techniques described in last section. There are three channels 
which contribute, T —)■ g*g*g* PeJ/i’, P —t g*g*l* —t PeJ/y, and T —)■ 7 * —)■ PeJ/i>- In 
the subsequent three subsections, we will calculate the SDCs of the three channels separately. 


A. T ^ g*g*g* ^ 7 , 

The typical Feynman diagram is illustrated in Fig. 1(a), which is drawn by using 
JaxoDraw [28]. According to the Feynman rules, we write down the amplitude of 
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FIG. 1: The typical Feynman diagrams for T —)> where Fig. 1(a) and Fig. 1(b) correspond 

to the channels T —>• g*g*g* —)• gcJ/'fp and T —)• g*g*^* gcJ/'4’^ respectively. 





FIG. 2: The typical Feynman diagrams for T —?■ 7 * —)• Fig. 2(a) corresponds to the 

tree-level diagram, and Fig. 2(b) and Fig. 2(c) correspond to the diagrams of the NLO radiative 
corrections from the initial state and final states, respectively. 


bb{po,po) —>■ g*g*g* —)■ cc{pi,pi) + cc{p 2 ,p 2 )- The corresponding spin states for the heavy 
quark pairs are readily obtained using (5), and the relativistic corrections are calculated 
by employing (7) and (8). In our calculation, we use the Mathematica package Fey- 
nArts [29] to generate the Feynman diagrams and Feynman amplitude, FeynCalc [30] 
as well as FeynCalcFormLink [31] to implement the algebras of Dirac trace and Lorentz 
indices contraction, and Fire [32, 33] and self-written programs [34] to make the tensor 
reduction. After some tedious work, we hnally get the amplitude, accordingly, the SDCs are 
obtained by removing the tensor and vectors in (1) and producing the matching factor 


1 


2N,2NAE{q,y2NAE{q2) 


2 ■ 


( 9 ) 


Since there are E{qi) in (9), we should further expand them by using (8). 

The SDCs are hnally obtained analytically. As we will see in the next section, the channel 
T —>■ g*g*g* PcJ/4’ dominates T decay into pc + hence, it is both necessary and 
valuable for us to present the explicit expressions for these SDCs in the Appendix. Our 
LO SDCs agree with those in [16]. The numerical predictions will be given in the following 
section. 
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B. T ^ ^ 


The Feynman diagram is illustrated in Fig. 1(b). For this channel, the J/'ip is produced 
through the fragmentation of a photon. At leading order in the channel can be splitted 
into T rjc'y* and 7 * —)■ J/ip. As discussed in Ref. [35], to reduce the theoretical uncertainty, 
the authors apply the VMD method [36] to calculate the fragmentation of the 7 * to J/ip 
in studying the process e+e“ ^ 7 * —)■ ?7c7* VcJ/4’- Similarly, here we adopt the same 

treatment as that in Ref. [35]. To this end, we first calculate the amplitude of the subprocess 
Vcl*, then multiply the amplitude with a factor 


gj/^Vina 

"'Ah 




where the effective coupling gj/^ reads [36] 


fi'j/h 




F[J/^ ^ e+e-; 


1/2 


( 10 ) 


( 11 ) 


where F[J/'^ —)■ e’''e“] signihes the decay rate of the process J/tp ^ e+e“. 

To get the SDCs, we apply the same techniques as in the last subsection. Finally, we 
obtain the analytic expressions for the amplitude, and immediately also for the SDCs. Unfor¬ 
tunately, the expressions of these SDCs are much cumbersome, therefore we do not attempt 
to present them in this paper, and merely yield the numerical predictions in the next section. 


T 
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gcJ/'P 


The Feynman diagrams are illustrated in Fig. 2. To the accuracy considered in the current 
work, the amplitude of this channel can be separated into two subprocesses, T —)■ 7 * and 
7 * —^ VcJ/ip- The NLO relativistic and radiative corrections to both subprocesses have been 
investigated. 

For T —)■ 7 *, the amplitude accurate up to O^asv"^) reads [37] ^ 


A(T ^ 7 *) = ^ 7 *) 


2a,C'p 


TT 


+ 


1 oigCp 

- + 

6 An 




. ( 12 ) 


where /ij represents the factorization scale. For 7 * —)■ the amplitude with the photon 

energy the same as that at B factories has been evaluated accurate up to Oi^aav"^) corrections 
in [39, 40]. In our process, since the virtual photon is produced through T annihilation, it 
is reasonable for us to take the invariant mass of 7 * to be the mass of T. After some effort, 
we are able to obtain the expression of the amplitude for the channel 7 * —)■ rjcJ accurate 
up to O^asy"^) corrections, and immediately yield the amplitude for T —)■ 7 * —)■ rjcJ/'tp with 
the help of (12). Our result for 7 * —)■ rjcJ/P^ with the same energy as B factories agrees with 
Refs. [39, 40], where the subprocess 7 * —)■ rjcJ/P^ has been deeply investigated in [39, 40]. 
For details, we refer the authors to Refs. [39, 40]. In the following section, we will present 
the numerical predictions. 


^ It can also be read from Ref. [38] by taking TOc = m{,. 
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IV. NUMERICAL PREDICTIONS AND DISCUSSIONS 

In this section, we select input parameters, and utilize the obtained analytic expressions 
to make numerical predictions for the SDCs as well as the decay rates. 

A. Choosing the input parameters 

To get the decay rates, we should calculate the squared amplitude and the phase space 
factor. The squared amplitude is directly derived by multiplying the amplitude in (1) with 
its complex conjugate. Sum over the polarizations for the final state J/'ip and average the 
spin states for the initial state T, we get 

^ = I ■ ^ 2 )" - PiPl] X l-4p. (13) 

spin 

where, to maintain the gauge invariance, one must take Pp = 4 E‘^{qi), P2 = 4i?^(g2) and 
Pi ■ P2 = 2 [T'^(go) — E‘^{qi) — E‘^{q2)]. On the other hand, the phase space factor does not 
affect the gauge invariance; therefore, we use the physical masses of the quarkonia there [37] 

(14) 

where A is defined via A (a, b, c) = + b‘^ + — 2ab — 2ac — 26c. 

We now specify our choices of the parameters in our computation. We take the heavy 
quark pole masses as rric = 1.4 GeV and mi, = 4.6 GeV. The masses of quarkonia are taken 


from Ref. [22]: 

mr{is) = 9.46030 GeV,mT( 2 S) = 10.02326 GeV,mT( 3 S) = 10.3552 GeV, (15a) 
= 2.9836 GeV,m^^( 2 s) = 3.6394 GeV, mj/^ = 3.096916 GeV. (15b) 

The matrix elements for J/ip are taken from Ref. [41], 

= 0.446 GeV^ = 0.223, (16) 

which are htted through the decay J/ip e^e~ accurate through the relative order and 

the matrix elements of rj^ are taken from Ref. [42], 

{0)n„iis) = 0.398 GeV^ {v^),,^iis) = 0.232, (17a) 

{0),.(25) = 0.202 GeV^ = 0.255, (17b) 

which are htted through the decay rjc^lS) —)■ 27 , ric{lS) —)■ LH, and ric{2S) LH. In (16) 


and (17), we have used the vacuum saturation approximation [1] to relate the production 
matrix elements with the decay ones. 

In addition, we should also choose the values of the matrix elements for T, which have 
been htted in Refs. [20, 43]. Here, we cite the values in Ref. [43], 


= 3.069 GeV^ 

{v‘^)r(is) — —0.009, 

(18a) 

(0)t( 2S) = 1.623 GeV^ 

{v‘^)r{2S) = 0.09, 

(18b) 

((!1)t( 3S) = 1.279 GeV^ 

{v^)ri3S) = 0.155, 

(18c) 
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which are fitted through the decay T e~^e~. 

In the channel T — )■ g*g*Y — )■ we also need the decay rate for J/ip e^e~, 

which is taken from the Ref. [22] as T[J/'ip e’''e“] = 5.55 keV. Since the hne structure 
constant runs slightly with the energy scale, we uniformly take a = ^. Finally, we choose 
the renormalization and factorization scales to he fj,R = Hf = nib, and the strong coupling 
constant asinih) = 0.22 [16]. 


B. Numerical predictions for the SDCs 

In this subsection, we present the SDCs numerically by using the parameters selected 
in the previous subsection. As mentioned, there are three channels which contribute to 
SDCs; hence, we add superscripts in the SDCs to make the distinction, i.e., where 
i = 1,2,3 corresponds to the channels T —)■ g*g*g* t g*g*l* VcJ/ip, and 

T —)■ 7 * —)■ rjcJ/'P’, respectively. Substituting the parameters into the analytic expressions of 


the SDCs, we obtain for T(IS') —)■ g*g*g* —)■ rjci)-S)J/'ip, 

4^^ = (-2.16 X 10"® - 8.90 X 10"^i)a^ GeV^ (19a) 

eg = (-5.64 - 1.83i), c^l = (-1.14 - 0.90i), = (-0.99 - 0.27f); (19b) 

for T(l^) ^ g*g*r ^ hc(l^) W, 

= (1.90 X 10"^ + 7.79 X GeV^ (20a) 

eg = (-3.34 - 0.62i), eg = (-0.73 + 0.31i), eg = (0.19 + 0.02i); (20b) 

and for T(IS') ^ 7* —)■ gciP^S)J/tp, 

eg = (4.14 X 10-® + 6.78 x lO’^a, - 7.08 x 10-^a,i)aa, GeV^ (21a) 

eg = (-0.17 - 0.05a,), eg = (0.52 - 1.60a, + 2.50a,i), (21b) 

eg = (0.35 - 0.89a,+ 1.76a,i). (21c) 


To see the dependence on a and a,, we retain these two constants in the SDCs explicitly. 
From (19-21), we notice some characteristics. First, T —)■ g*g*g* —)■ gcJ is the dominant 
channel for the SDCs and therefore also for the decay rates. Second, the SDCs of the 
three channels are of quite different phase, in addition, the phases of various relativistic 
corrections are also distinguishing. The last is not the least, the relativistic corrections from 
T are extraordinarily large, so one may expect that this corrections will also contribute a 
lot to the decay rates even though the corresponding matrix element is rather small. 


C. Numerical predictions for the decay rates 

Combining the squared amplitude (13) with the phase space factor (14), we are able to 
evaluate the numerical predictions for the decay rates of various processes. All the decay 
rates are listed in Table. I. 

In the table, we give the LO decay rates Flo as well as the contributions from the 
relativistic corrections. To see it clearly, we separately list the relativistic corrections from 



TABLE I: The LO and NLO relativistic corrections to decay rates of T —)■ rjcJ/i^- Tlq represents 

(i) 

the decay rate at leading order in velocity expansion; Tj^j^o i = 0,1,2 represent the relativistic 
corrections from T, rjc, and J/V', respectively; Tnlo represents the total relativistic corrections. 


Processes 

Plo (eV) 

(eV) 

^nIo (eV) 

^'nIo (eV) 

^NLO 

Tt.o 

T(1F) ^ pc(15) J/V- 

0.41 

0.03 

-0.28 

-0.22 

-115% 

T(1F) ^ Pc(25) J/V- 

0.23 

0.02 

-0.17 

-0.13 

-122% 

T{2S) ^ gc{lS)J/^ 

0.22 

-0.18 

-0.15 

-0.12 

-204% 

T{2S) ^ Pc(25) J/V- 

0.13 

-0.11 

-0.10 

-0.07 

-211% 

T(3F) ^ gc{lS)J/^ 

0.18 

-0.25 

-0.12 

-0.10 

-263% 

T(3F) ^ Pc(25) J/V’ 

0.10 

-0.15 

-0.08 

-0.06 

-270% 


different sources; i.e., with i = 0,1,2 originates from T, rjc, and J/'^, respectively. 

From the table, we hnd that almost all the relativistic corrections are both signihcant and 
negative, except for F^lo Tf'(lS') decay due to a tiny and negative matrix element given in 
(18). Both the relativistic corrections from r]c and J/ijj are negative, which is quite different 
with that in rjc associated with J/ijj production at B factories [44]. The discrepancy may be 
accounted for by distinct Feynman diagram structure. 

Another distinctive feature in the table is that the relativistic corrections from T(2S', ?>S) 
are extraordinarily large, which even outnumber the corrections from charmonia, though the 
NLO NRQCD matrix elements are smaller than these for charmonia. We may learn a lesson 
from this kind of process and the example T —)■ ccgg in Refs. [20, 21]. Despite a rather 
small NRQCD matrix element (n^)T, it is not always reasonable to discard the relativistic 
corrections from T, which actually have not been considered carefully in most references. 
Moreover, it is also important and valuable to further determine the corresponding NRQCD 
matrix elements for T through independent approaches. 

As mentioned in the introduction, T —)■ rjcJ is of a similar Feynman diagram structure 
as pn. We may expect there are also a relative large relativistic corrections for 'ijj in 
tfj —)■ pn. Nevertheless, due to both large (even huge) and negative corrections, the “p — vr” 
puzzle is unlikely to be explained purely by the relativistic corrections, or at least by the 
NLO corrections. 

By combining all the corrections, we show the ratios of the decay rates of the relativistic 
corrections to that at leading order in the Table. I. Our results indicate that the NLO 
relativistic corrections exceed the LO contributions for all the processes, and what is worse, 
the corrections are negative, which renders the total decay rates negative ^ and therefore 
unpredictable at next-to-leading order. A resummation for the large relativistic corrections 
may be needed [37, 38], but is, however, out of the scope of this work. 

In summary, we study the NLO relativistic corrections to T exclusive decay into double S- 
wave charmonia. We calculate the contributions from three channels T —)■ g*g*g* —)■ Pc+J/'4’i 
T —)■ g*g*p* —)■ Pc + J/i^i and T —)■ 7 * —)■ pc + J/i^- Our results indicate, in spite of the 
small (n^)T, the relativistic corrections from T are extraordinarily large, which suggests the 
relativistic corrections from T are not always negligible. Moreover, we hnd that all the 


^ Since decay rates are expanded in powers of v and as, and truncated at relative 0(v^), we obtain negative 
values. 
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relativistic corrections to the decay rates are both signihcant and negative, which indicates 
a poor convergence in velocity expansion for this kind of process. 
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Appendix A: Analytic expressions of SDCs cO for the channel T 

For convenience, we rewrite the SDCs as 


9*9*9* 


VcJ/'ip 


Ji) _ ^ 
° Qml 


a. 


a" 


Cq C 2 1 *"0 ^ 2,2 ~ 


6ml 


* 2 , Cq C 2 3 — 


at 


Qml 


■^3: 


(Al) 


and dehne 


LI 

L2 

L4 

L8 

II 

Al 


Re 

Re 


(_1 + ,V3/3) In iUzil + (1 + ,^131 In *('3 + 1) 


V 3 / 3 - 

j. {-i + \/ 3 )l 3 {-i + \/ 3 )l 3 

LI2-;=- LI2 


Vsp + i 




, L3 = Li 


VS/S + i\ ’ 
1-/3 


T- /3-1 T. T- 2/3 T. T- T. T- 

L12 ^ 1 , L 5 — L12-5——, L6 — L12-5——, L 7 — L12- 


/3 + 1 


^/3 + r 


/3-1 


' 2 - 6r ’ 


u/4^. L9 = LIO = U.'UA, Lll = Li. 


'4(1 - 3r) 

In 2,12 = ln(l — /3), 13 = In 
tan“^(\/3/3). 


/3 + 1 
1-/3 


1 + /3 


6r — 2 
, 14 = ln(l — 3r), 15 = In 


12r-4 ’ 
4(1 - 3r) 

(/3 + l)^’ 


(A2) 


where r 



, (3 = 1/1 — 4r, and Li 2 (a;) represents the dilogarithm. 


10 



Utilizing the definitions in (A2), we present sq, si, S 2 , and S 3 as follows 


-Sl 


-So = 


207r(2r - l)r3(Il^ + 12^ + 2L3 + L4) 57r(10r - 3)r3l32 


27/33 27/33 

57r(2r — l)r[8r^I3 + /3(2r — 1)]I2 57r(2r — l)r(/3 + 8r^I2 — 8r^I3 — 2/3r)Il 


57rr 


+ 


27/33 27/33 

/3 — 8(214 + m + 4)r3 + 4(/3 + ivr + 4)r^ — 2(2/3 + 3)r + 1 


13 


407rr"^ / 27r 


54/33 

-A1 + 2L2 + 2L5 + L6 + L7 - L8 - L9 - LIO + Ill 


27/33 V 3 

57rr^ 3/33 + 2(2\/3 — 3i)7r/3r + n‘^r{2r — 1) 


8 I /33 


(A3) 


407rr3Ll 407rr^(I4 — 15) 107r(104r^ — 8r — 9)r3 27r 


+ 


27(3r - l)/35 27/34 

-L8-L9-L10 + L11) + 


8 I /35 

107r(304r2 - 206r + 37)r3 


A1 + 2L2 + 2L5 + L6 + L7 


8 I /35 


(11^ + 12^ + 2L3 + L4) 


57r(1120r2 - 634r + 93)r3l3^ 


162/3^ 


57r 


972(12r2 - 7r + 1)2/3 


27r2(l - 3r)2(304r2 - 206r + 37)r^ + Att 


8(-405i + 18iy3)r^ 


+ (3213i - 1532y3)r2 + 6(-177i + 88y3)r - 60^3 + 117i 


(3r^ + (-14976r® + 12888r® 


-1162r^ - 2593r3 + 1286r2 - 249r + 18)/3 \ + 


57r 


24(304r2 - 206r + 37) 


486/33r 

r^(1113 - 1112 - 1213) + (480r® - 44r^ - 200r3 + 235r2 - 84r + 9)(II - I2)/3 

57ri3 r <93 r 

+48^1 (i2;r^WTWTT) 

-48i(7/3 + 33))r® + 4(14598/3 + 3i7r(1839/3 + 5585) + 34606)r^ + (-20974/3 
-6i7r(1873/3 + 3585) - 25754)r® + (-5731/3 + 24i7r(107/3 + 144) - 18253)r^ 
+2(4749/3 - llli7r(/3 + 1) + 8515)r^ - 2(2353/3 + 3431)r=^ + (1225/3 + 1537)r2 

-3(55/3 + 61)r + 9(/3 + 1) - 12r^(104r2 - 8r - 9)I4l, 


(A4) 


11 



S2 


107r(13r — l)r^ / 2 ti 


8 I /35 


-A1 + 2L2 + 2L5 + L6 + L7 - L8 - L9 - LIO + Lll 


207r(40r3 - 15r + + 12^ + 2L3 + L4) 57r(120r3 - 32r + 5)r3l3^ 


Stt 


8I/35 8I/35 

48(40r^ - 15r + 2)r3(IlI3 - 1112 - 1213) + /3(288r® - 352r^ 

207rr^(I4 - 15) 207rr^Ll 


486/^5 

+224r2-48r + 3)(Il -12) 
57rI3 


232r^ 


27/3^ 


- 12(13r - l)r'^I4 - 


/S' 


27/35(3r - 1) 

— 3(/3 + 1) 


486/35 [ ' ' (/3 +l)(12r2 - 7r + 1)2 

-384i(457r - 236i)r® + 96(290/3 + 15i7r(3/3 + 11) + 1402)r® 

+16(-2523/3 - 1537r(12/3 - 7) - 5003)r^ + 4(3842/3 - 3i7r(95/3 + 527) + 2110)r® 
+4(1742/3 + 324i7r(/3 + 2) + 5075)r® + (-8805/3 - 12i7r(27/3 + 35) - 14809)r^ 

+ (3583/3 + 24i7r(/3 + 1) + 4913)r3 - 2(367/3 + 439)r2 + (75/3 + 81)r I 


hnr 


3671^(1 - 3r)2(40r' - 15r + 2)r^ 


8748/3(12r2 - 7r + 1)2 [ 

+477/3(216(70^3 - 937)r^ - 4(4036^3 - 440l7)r3 + (6468^3 - 593l7)r2 
-18(68^3 - 597)r + lOOVs - - 9/3(10560r‘^ - 9568r^ + 648r^ + 1983r3 


-824r2 + 123r - 6) , 


(A5) 


12 



and 


lOOTrr^Ll 1007rr^(I4 — 15) 107r(72r^ + 31r — l)r^/27r 

5*^ —-;-^—r — -:--::-I —A1 “h 21^2 2LjO IjO 

27(3r-l)/35 27f3^ 81f3^ 

+L7 - L8 - L9 - LIO + Lll^ 

, 107r(400r3 - 24r2 - 140r + 19)r3(Il2 + 12^ + 2L3 + L4) 

^ 81^5 

57r(1200r3 + 72r^ - 358r + 55)r3l3^ 


Stt 


187r2(l-3r)^ 


162/35 8748(12r2 - 7r + 1)2/3 

(400r^ - 24r2 - 140r + 19)r^ + 87r(540(-933 + 70V3)r^ + (435243 - 40036V3)r^ 
+3(-47943 + 529iy3)r2 - 90(-283 + 33y3)r + 241^3 - 234i)/3r^ - 9(52800r^ 

-57632r® + 22872r® - 7217r^ + 4017r3 - 1565r2 + 279r - 18)/3 


5% 


24(400r^ - 24r2 - 140r + 19)r^(IlI3 - 1213 - 1112) + (1440r® - 1904r® 


486/35r 

-4r^ - 100r=^ + 307r2 - 99r + 9) (II - I2)/3 
57rI3 f /3^ 


10 


486/35r \ (12r2 - 7r + If {(5 + 1) 

-963(371(75/3 + 293) - 23(725/3 + 3421))r9 + 16(12525/3 + 337r(327/3 - 31) + 25513)r® 
+83(53 771/3 + 127973/3 + 349871 + 163493)r^ + 2(11339/3 - 337i(987/3 + 2003) + 15363)r® 
+23(76271/3 + 52403/3 + 99071 + 126313)r^ + (11041/3 - 11437i(/3 + 1) + 20095)r^ 

-4(1411/3 + 2046)r3 + 2(716/3 + 887)r2 - 18(10/3 + ll)r + 9(/3 + 1) 


-12r^(72r2 + 31r - 1)14 


(A6) 
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